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The critical effect of disorder on the two-dimensional (2D) surface superconductor Si(lll)- 
(^/7x^/3)-In is clarified by comparing two regions with different degrees of disorder. Low- 
temperature scanning tunneling microscopy measurements reveal that superconductivity is retained 
in the less disordered region, judging from the characteristic differential conductance (dl/dV) spec¬ 
tra and from the formation of vortices under magnetic fields. In striking contrast, the absence of 
those features in the highly disordered region shows that superconductivity is strongly suppressed 
there. Analysis of observed zero-bias anomalies in dl/dV spectra allows us to estimate the reduction 
in the transition temperature Tc, which explains the fate of superconductivity in each region. 


Recently, the emergence of superconductivity has been 
established for several kinds of silicon surface reconstruc¬ 
tions with metal adatoms [1-5]. They are particularly 
interesting because of their well-defined unique atomic 
structures and accessibility through standard surface sci¬ 
ence techniques such as scanning tunneling microscopy 
(STM) [1, 6, 7]. One of the important features of these 
systems is a high sensitivity to the presence of surface 
defects, which originates from the atomic-scale thick¬ 
ness of the conducting layer. For example, the influ¬ 
ence of isolated surface atomic steps on supercurrents 
was studied by electron transport measurements [2] and 
with a low-temperature (LT) STM under magnetic fields 
[6, 7]. Detailed analyses in these experiments have re¬ 
vealed that the atomic steps can play the role of Joseph- 
son junctions. In contrast, a high concentration of ran¬ 
domly distributed atomic steps and point defects should 
be regarded as crystalline disorder. In this case, the 
disorder-induced suppression of superconductivity is ex¬ 
pected. Although Anderson’s theory claims that super¬ 
conductivity is insensitive to disorder under time-reversal 
symmetry [8], electron localization inherent to the disor¬ 
dered two-dimensional (2D) system [9] could lead to a 
strong suppression of superconductivity [10-12]. More¬ 
over, STM studies on conventional thin films of super¬ 
conductor have shown that the presence of strong disor¬ 
der makes the superconducting state spatially inhomo¬ 
geneous [13-15]. However, such effects still need to be 
investigated for surface 2D superconductors. 

In the present Rapid Communication, we successfully 
clarify the critical influence of disorder on the supercon¬ 
ductivity in the Si(lll)-(V^xv/3)-In surface reconstruc¬ 
tion [referred to as xv/3)-In] using LT STM. Topo¬ 
graphic STM observations show two surface regions with 
different degrees of disorder. In the less disordered re¬ 
gion, superconductivity is found to survive judging from 
the characteristic differential conductance (dIjdV) spec¬ 
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tra and from the formation of vortices under magnetic 
fields. In striking contrast, these features are absent 
in the highly disordered region, showing strong suppres¬ 
sion of superconductivity. Analysis of observed zero-bias 
anomalies in dl/dV spectra allows us to estimate the 
degree of disorder and the reduction in the transition 
temperature Tc. The result satisfactorily explains the 
survival and disappearance of superconductivity found 
in each region. 

The (7? x\/3)-In surface was chosen as an archetypal 
superconductor made of a surface reconstruction [1-5, 7]. 
Although the existence of two (\/7 x\/3)-In phases was 
reported [16], the surface atomic structure observed here 
is only of one kind and is identical to that in our previous 
studies [2-4, 7] [17]. In the following, the physical param¬ 
eters needed for analysis will be taken from the angle- 
resolved photoemission spectroscopy (ARPES) study on 
a (77 x\/3)-In sample prepared in the same way as ours 
[18], with an effective mass of m* = 1.09me and a Fermi 
energy of Tp = 6.9 eV. These give the Fermi wave num¬ 
ber = 1.41 and the Fermi velocity 'i;F = 1.49 x 10^ 
m/s. 

The experiments were performed using a LT STM sys¬ 
tem. First a Si(lll) substrate was flashed several times 
to prepare a clean surface with the 7x7 reconstruction. 
Indium was thermally deposited onto the surface at room 
temperature and annealed to obtain the ^/l x^/3 phase. 
Here we employed an annealing condition of ^400 °C for 
5 min, which was lower than that of ^570 °C for 3 s 
used in the previous study [7]. This decrease in anneal¬ 
ing temperature resulted in a substantial increase in the 
density of surface defects. The presence of the ^/7 x^/3 
phase was confirmed by reflection high-energy electron 
diffraction (RHEED) patterns. The sample was then 
transferred to the LT STM stage where it was cooled 
down to ^0.5 K, which is sufficiently lower than Tc ^ 3 
K [1-5]. Differential conductance dl/dV (/: tunneling 
current; V: sample bias voltage) was measured by us¬ 
ing a lock-in amplifier with a small ac bias modulation. 
Zero-bias conductance (ZBC), i.e., dl/dV at F = 0, was 
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FIG. 1. (Color online) (a) STM image of a 500 nm x 500 
nm region on the surface (setpoint: 10 pA at +50 mV), (b) 
STM images of a 20 nm x 20 nm region on a flat region (set- 
point: 50 pA at +500 mV). Note that the void defects in this 
image appear as white dots in (a) and in Fig. 2(b) because of 
insufficient spatial resolutions in the latter images, (c) STM 
images of a 20 nm x 20 nm region on a rough region (set- 
point: 50 pA at +500 mV). The arrows indicate the holes 
where In films did not grow. The parallelograms in (b) and 
(c) represent y/7 Xy/S unit cells, (d), (e) Fourier transforms 
of the STM images (b) and (c), respectively. The circles indi¬ 
cate spots corresponding to the y/7 X\/3 periodicity, (f) Line 
prohles measured along the solid lines in the STM image (a). 


measured while the feedback was off after the tip-sample 
separation was fixed. Magnetic fields were applied in the 
direction perpendicular to the sample surface using a su¬ 
perconducting solenoid magnet. 

First, we characterize the morphology of the sample 
surface. Figure 1(a) shows a representative STM im¬ 
age in which the fiat central area and the surrounding 
defective region are visible. In the following, they are 
referred to as “fiat” and “rough” regions, respectively. 
An atomic-resolution image of the fiat region reveals the 
presence of a well-ordered x^/3 reconstruction, the 
unit cell of which is indicated as the parallelogram [Fig. 
1(b)]. The image also includes a low density of small void 
defects; the average spacing between the voids was found 
to be nm from a larger STM image. Accordingly, the 
Fourier transform (FT) of the image exhibits clear spots 
corresponding to the x^/3 periodicity [Fig. 1(d)]. In 
contrast, the rough region is characterized by a high den¬ 
sity of defects as evidenced by an atomic-resolution im¬ 


age with the same scale [Fig. 1(c)]. Nevertheless, the 
\/7 x\/3 reconstruction is locally preserved within small 
domains 5-10 nm in size (see the parallelograms). The 
presence of the \/7 x\/3 periodicity is also confirmed by 
its FT pattern, in which the corresponding weak spots 
are visible [Fig. 1(e)]. The line profiles taken across 
the domains reveal that they differ in height by 0.3 nm 
[Fig. 1(f)], indicating that the domains are separated by 
atomic steps of silicon. There are often holes along the 
step edges, as shown by the arrows in Fig. 1(c). Such 
structures are commonly observed along the atomic steps 
of the (^/7 x\/3)-In surface [7]. To summarize, both fiat 
and rough regions are \/7 x\/3 surfaces, but with very dif¬ 
ferent degrees of crystalline disorder. This atomic-scale 
characterization of disorder became possible because the 
(\/7 Xv/3)-In superconductor has a clean surface that is 
suitable for STM. As shown in the following, the ob¬ 
served morphological differences have a profound effect 
on superconductivity. 

Figure 2(a) depicts a series of dl/dV spectra taken 
across a boundary between the flat and rough regions. 
The spectral locations are indicated in the STM image of 
Fig. 2(b). The dl/dV spectra taken in the fiat region {A- 
B) are characteristic of the superconducting state; each 
spectrum exhibits an energy gap at V = 0 and two co¬ 
herence peaks at the gap edges. Despite the presence 
of void defects, the magnitude of the energy gap was 
found to be almost constant within the flat region. Fig¬ 
ure 2(c) shows a representative spectrum in the flat re¬ 
gion measured with a higher energy resolution. Fitting 
the spectrum to the Dynes formula [19] yields an energy 
gap A = 0.45 meV, a broadening parameter F = 0.08 
meV, and an electron temperature T = 1.28 K. The ob¬ 
tained A is ^20% smaller than the value of 0.57 meV 
reported previously [1]. The reason for this discrepancy 
will be discussed later. In contrast, the dl/dV spectra 
taken in the rough region exhibit only a much broader 
dip structure around V = 0 [see C-D in Fig. 2(a)]. The 
absence of the coherence peaks in these spectra indicates 
the disappearance of superconductivity in the rough re¬ 
gion. 

The above findings were corroborated by observing 
ZBC images of the same area in different magnetic fields 
[Figs. 3(a)-3(d)]. At the lowest field of 5ext = 0.08 T, 
two bright spots corresponding to high ZBC were ob¬ 
served in the flat region [Fig. 3(a)] and they increased in 
number and overlapped for 5ext =0.12 and 0.16 T [Figs. 
3(b) and 3(c)]. Thus the bright spots are attributed to 
vortex cores [20-24] and superconductivity is evidenced 
in the flat region. At 5ext = 0.24 T, ZBC was saturated 
over the flat region, meaning that the superconductiv¬ 
ity was completely suppressed [Fig. 3(d)]. On the con¬ 
trary, ZBC was not affected by these magnetic fields in 
the rough region. This is clearly visible in Fig. 3(e), 
where the averaged ZBC is plotted as a function of dis¬ 
tance d from the boundary between the two regions for 
each magnetic field. While ZBC increases with increas¬ 
ing magnetic field in the flat region {d > 0), reflecting 
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FIG. 2. (Color online) (a) Series of tunneling spectra mea¬ 
sured in a zero magnetic field at 15 points indicated in the 
STM image in (b). The curves are shifted from each other 
by 4 nS for clarity (setpoint: 400 pA at +20 mV; bias mod¬ 
ulation: 200 /iVrms at 610 Hz), (b) STM image of a 50 nm 
X 200 nm region including a boundary between a flat region 
and a rough region (setpoint: 10 pA at +50 mV). The white 
dots in the flat region correspond to the void defects in Fig. 
1(b). (c) Representative tunneling spectrum of the flat region. 
The dashed curve is the result of fitting by the Dynes formula 
(setpoint: 400 pA at +20 mV; bias modulation: 50 /xVrms at 
610 Hz,), (d) Semilogarithmic plot of the tunneling spectra 
at points A and D. The solid (dotted) curves represent the 
data on the positive (negative) biases. The black solid lines 
are the results of fitting by Eq. (1). 





FIG. 3. (Color online) (a)-(d) ZBC images in magnetic 
fields of Hext = 0.08, 0.12, 0.18, and 0.24 T. The field of view 
is the same as that of Fig. 1(a). Vortex cores are imaged 
as bright spots in the flat region (setpoint: 200 pA at +20 
mV; bias modulation: 200 /iVrms at 610 Hz), (e) Average 
ZBC plotted as a function of distance d from the boundary of 
the flat region. To plot these curves, ZBC was averaged over 
all the points with a given value of d, where the sign of d is 
defined to be positive (negative) for the flat (rough) region, 
(f) Line profile of the ZBC image along the line P-Q across a 
vortex core in (a). 


the suppression of superconductivity in the vortex cores, 
ZBC is almost constant in the rough region (d < 0). This 
strongly indicates that the superconductivity was already 
suppressed in the rough region and that the dip structure 
observed for C-D in Fig. 2(a) cannot be attributed to 
superconductivity. 

Considering the high degree of disorder in the rough 
region, this type of dip structure is ascribed to the zero- 
bias anomaly (ZBA) originating from electron-electron 
interactions enhanced by electron scatterings [25-27]. In 
2D systems, the ZB A is known to exhibit a logarithmic 
bias dependence [15, 28-35]. According to perturbation 
theory calculations [26, 27], the correction to the density 
of states 6p to the unperturbed value p = ra "due 


to this effect has the following form, 


P 


= 1 - 


dTr/cp/ 


In 




( 1 ) 


with 5 = max{|eV|, k^T}. Here, I is the elastic mean free 
path, r = 1/vy the elastic scattering time, D = {1/2)v-fI 
the 2D diffusion constant, and k = e^pj{2eo) the 2D 
inverse screening length. In the energy range of 1-20 
meV, the first logarithmic factor of the right hand side is 
roughly constant, and the ZBA has a ln(V) bias depen¬ 
dence coming mainly from the second logarithmic fac¬ 
tor [36]. Figure 2(d) replots the dl/dV spectra taken 
well inside the rough region (point D) and the flat region 
(point A). The spectrum at D exhibits a clear ln(V) bias 
dependence above 1 mV, which supports our interpreta- 
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TABLE L Parameters determined from fitting analyses for 
the spectra shown in Fig. 2(a). The errors come from the 
standard deviations of I values determined at various posi¬ 
tions. 



A (meV) 

1 (nm) kpl 

Tc (K) 

? (nm) 

Flat region 

0.45 

2.07 ±0.56 29 ±8 

1 4+0-4 

4- + -0.6 

38 ±5 

Rough region 

NA 

1.05 ± 0.12 15 ±2 0.19lo;l3 

NA 


tion of the dip structure as a ZB A. Assuming k-p = 1.41 
and Up = 1.49 x 10^ m/s as mentioned earlier [18], 
/ = 1.05 ± 0.12 nm and k-pl = 15 ± 2 are obtained by fit¬ 
ting Eq. (1) to the spectrum (fitting range: 2 mV <V< 
20 mV). Here, k^l is the dimensionless measure of the de¬ 
gree of disorder. Thus the rough region is metallic in the 
sense that k-pl ^ 1. This is consistent with the absence of 
a Coulomb gap in the spectra, which should be observed 
for an insulating state with k^l 1 [32, 33, 37, 38]. We 
note that the dl/dV taken at A in the flat region also 
shows a ln(V) dependence above the superconducting co¬ 
herence peaks, although its dependence is much weaker. 
This means that a weak disorder-induced ZBA coexists 
with superconductivity in the flat region. A similar fit¬ 
ting analysis gives I = 2.07 ± 0.56 nm and k-pl = 29 ± 8 , 
which are larger than those for D. The obtained values 
of I and k^l are summarized in Table I, together with 
other parameters. 

The strong suppression of superconductivity in the 
rough region can be explained in terms of the combined 
effects of electron-electron interactions and disorder [39- 
41]. According to the theory, the decrease in Tc is given 
by [41] 


A = exp(-l/ 7 ) 


X 



Vij2 \ 


Vi/2 

'y-t/i 


j-| 1/V2t 


( 2 ) 


where Tco is Tc in the absence of disorder, 7 = 

1/ln(/cBTcoT/?i), and t = l/(7r/cF/). Using k^l obtained 
above and Tco = 3 K, Tc for the rough region is estimated 
to be 0.19tQ }3 K. This value is sufficiently lower than 


the experimental temperature of T ^ 0.5 K. Thus the 
absence of superconductivity in the rough region is ra¬ 
tionalized. In contrast, Tc for the flat region is estimated 
to be IA'^q'q K from a similar argument, consistent with 
the presence of superconductivity in that region. The 
finding is also in line with the observed suppression of 
the superconducting energy gap A by ^20%, since Tc 
and A are proportional in the Bardeen-Cooper-Schrieffer 
(BCS) theory. The discrepancy between the estimated 
reductions in Tc and A may reflect the limitation of the 
perturbation approach used in the ZBA theory [25-27] 
and/or the deviation of the surface states from the nearly 
free electron model. 

Lastly, we estimate the influence of disorder on the 
coherence length ^ in the flat region. In a clean super¬ 
conductor, ^ at T = 0 is approximately equal to the 
BCS value defined by ^0 = hv^/iTrA). Here, ^0 = 694 
nm is obtained from A = 0.45 meV in the flat region 
and Up = 1.49 x 10^ m/s. In the presence of disorder- 
induced scatterings, ^ is given by ^ ^ Adopting 

I = 2.07 ± 0.56 in the flat region leads to ^ = 38 ± 5 nm. 
^ is independently determined to be 36 nm from the half 
width at half maximum of the ZBC profile measured over 
a vortex core [Fig. 3(f)]. The fact that these two values 
of ^ are nearly equal demonstrates the consistency of our 
analysis. 

In conclusion, we have revealed that superconductiv¬ 
ity was retained in the flat region of the \/7 X\/3 surfaces, 
judging from the presence of an energy gap and coher¬ 
ence peaks in the dl/dV spectra and from the formation 
of vortices under magnetic fields. In striking contrast, the 
absence of those features in the rough region showed that 
superconductivity was strongly suppressed there. Analy¬ 
sis of the zero-bias anomalies in dl/dV spectra allowed us 
to estimate the reduction in transition temperature Tc, 
which explained the fate of superconductivity in each re¬ 
gion. The present finding demonstrates the applicability 
of the existing theories to this class of surface 2D super¬ 
conductors, and lays the groundwork for future studies 
on the influence of surface defects. 
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